Universality in Complex Wishart ensembles: The 2 cut 

case 



M. Y. Mo 

oo 
O 

o 

Abstract 

' We considered N x N Wishart ensembles in the class Wq (£tv , M) (complex 

Wishart matrices with M degrees of freedom and covariance matrix Sjy) such that 
psj ', No eigenvalues of is 1 and Ni = N — Nq of them are a. We studied the limit as 

M, N, N and Ni all go to infinity such that jfe — > c, ^ — >■ /5 and < c, /3 < 1. In 
this case, the limiting eigenvalue density can either be supported on 1 or 2 disjoint 
Q H intervals in K + . We obtained a necessary and sufficient condition on the parameters 

a, (3 and c such that the limiting distribution is supported on 2 disjoint intervals and 
have computed the eigenvalue density in the limit. Furthermore, by using Riemann- 
Hilbert analysis similar to the one in [TT] (See also [27]), we have shown that under 
proper rescaling of the eigenvalues, the limiting correlation kernel is given by the sine 
kernel and the Airy kernel in the bulk and the edge of the spectrum respectively. As 
a consequence, the behavior of the largest eigenvalue in this model is described by 
<0 the Tracy- Widom distribution. 
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Let X be an M x N (assuming M > N) matrix with i.i.d. complex Gaussian entries whose 
real and imaginary parts have variance | and zero mean. Let S^r be an x N positive 
definite Hermitian matrix with eigenvalues ai,...,a^ (not necessarily distinct). In this 
paper, we will consider the case where E^r has only 2 distinct eigenvalues, 1 and a such 
that A^i of its eigenvalues are a and N — Ni of them are 1. We will assume that — > c and 



M 

j± — > f3 as N, M — > oo and that < c, (3 < 1. To be precise, we will assume the following 



cM- N = n = 0(1), N(3-N 1 = t 2 = 0(1), M,Ar,A^!^oo. (1.1) 

i i 
Let be any Hermitian square root of Sjy. Then the columns of the matrix XS^- are 

random vectors with variances |-v/%- Let the matrix be the following 

B N = ^IX^X^i, (1.2) 
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Then is the sample covariance matrix of the columns of XS^, while is the covariance 
matrix. In particular, Bn is a complex Wishart matrix in the class Wc (£jv, M). 

The sample covariance matrix is a fundamental tool in the studies of multivariate 
statistics and its distribution is already known to Wishart at around 1928 (See e.g. [2"S"] ) 



V{B N ) = ie"^" 1 ^ (det S) 



M—N 



M > N, 
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for some normalization constant C. 

Let yi > . . . > yN > be the eigenvalues of the sample covariance matrix B^ and let 
A, A be diagonal matrices such that 

A = diag(ai, . . . , a N ), A = diag(yi, ...,y N ), 

B N = U B AU B \ Sjv = U A AUJ\ U A , U B e C/(iV). 



Then the eigenvalue density of B^ is given by 



i<j 3=1 JQ 



-Mtr(A~ 1 QAQ- 1 ) 



dQ, 



;i-4) 



'QeU(TV) 



where Q = Ub € U(N) is unitary and c/Q is the normalized Haar measure on U(N). 
The eigenvalue density of -Bjv can be simplified by using the Harish-Chandra [23] (or 
Itzykson-Zuber [21]) formula 



/ 

JQ. 



N 



det 



-Mtr(A 



&U{N) 



- lQAQ -^dQ = l[T(j) 



-Ma: 



Vi 



J l<i,j<N 



3=1 



det 



-Ma, 



J l<i,j<N 
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where in the case of A or A having multiple eigenvalues, the above should be interpreted 
by using the L'Hopital rule. 

By using ( 11. 5ft in ( 11.41) . we see that in the case where N\ eigenvalues of is a and 
N — Ni is 1, the joint probability density function (j.p.d.f) for the eigenvalues of B^ is 
given by 



N 



Hv) 



J M.N 



y dj-i e _ MaJ i yk 



i<3 



3=1 



l<j,k,<N 
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where Zm,n is a normalization constant and dj , dj are given by 



d?=j, a j = l, l<j<N-N u 

d? = j-N + N 1 , aj =a, N — Ni < j < N. 

In this paper we will study the asymptotic limit of the Wishart distribution with — > c 

A' 



and 4jr — > /3 as M, iV — > oo in such a way that < j3, c < 1. In this case, the empirical 
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distribution function (e.d.f) Fn of the eigenvalues will converge weakly to a nonrandom 
p.d.f. F, which will be supported on either 1 or 2 intervals in R + . By applying the results 
of [H] to our case, we can extract properties of the measure F from the solutions of an 
algebraic equation (See Section [3 for details) 

za^ + (A 2 z + B 2 )i 2 + (z + + 1 = 0, 
A 2 = (l + a), B 2 = a(l-c), (1.7) 
B 1 = l-c(l-/3) +a(l-c0). 

Our first main result involves a necessary and sufficient condition for F to be supported 
on 2 intervals and an explicit formula for the distribution F (See Theorem [5]). 

Theorem 1. Let A be the discriminant of the quartic polynomial 

a 2 (l - c)e + 2(a 2 (l - c(3) + a(l - c(l - /3))£ 3 

+ (1 - c(l - /3) + a 2 (l - c/3) + 4a)£ 2 + 2(1 + a)£ + 1 = 0. 1 ' j 

i/ien t/ie support of F consists of 2 disjoint intervals if and only if A > 0. 

When A > 0, we also have the following expression for the density function of F (See 
Theorem E]) . 

Theorem 2. Let A be the discriminant of the quartic polynomial ( fl.ffj) . Suppose A > 
and Zet 71 < . . . < 74 6e i/ie ^ rea/ roots to (T772J). Let Ai < . . . < A4 be the following 

1 l-(3 ap , 1 
A fc = + C T~ + c 7^ > fc = l,...,4. 

7fc 1 + 7fc 1 + a 7fc 

Then the p.d.f F is supported on [Ai, A 2 ] U [A3, A 4 ] with the following density dF(z) = p(z)dz 
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2^ 



r ( z ) + \l-ztoKt D *( z ) \ ( r ( z ) - V -27^^3(2:) 



1.9) 



where D%(z) and r(z) are given by 



D 3 (z) = (l-a) 2 H(z~X 3 ), 

i=i 

1 / 2.Bf , / 9B 1 B 2 
r(z) = — [ ^z" 3 + ' 



+ 



271 a 3 

9A 2 2A 3 
a 2 a 3 



27 6A 2 5 2 . , 
a a° 
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The constants A\, B\ and B 2 in the above equation are defined by 



A 2 = (l + a), B 2 = a(l-c), B x = 1 - c(l - (5) + a(l - c(5). 



The cube root in $3.26}) is chosen such that y A e M for A G IR and the square root is 
chosen such that \fA > for A > 0. 

Remark 1. Since D 3 (z) can be written as 

lUZ) -27 \r{z)f- 4 {p{z))\ (1.10) 



a 4 z A 



for some polynomial p(z) in z~ l . We see that ifr(z) vanishes at any of the Xk, then D 3 (z) 
will have at least a double root at these points, hence r(Xk) 7^ 0. From this and $1.91) , we 
see that the density p(z) vanishes like a square root at the points X k , k = 1, . . . ,4. 

p{z) = ^\z-X k fi + 0((z- A fc )) , z -> X k . (1.11) 

7T 

Having obtained the global statistics of eigenvalues, we will continue to answer questions 
about the local statistics for the eigenvalues. In particular, we will show that, under 
suitable scaling of the eigenvalues, the behavior of the largest eigenvalue is given by the 
Tracy- Widom distribution. 

A result by Baik, Ben-Arous and Peche [7] shows that the correlation functions of the 
eigenvalues can be expressed in terms of a Fredholm determinant with kernel Km,n(x,v)- 
In [13] and [15], the authors have expressed this kernel in terms of multiple orthogonal poly- 
nomials (See Section [2] for details) and have shown that the m-point correlation function 
for the Wishart distribution (jl.3p is given by 

TZ^ N \ yi , ...,y m )= det (^(j/i,^))^ (1.12) 
where !Z^f ,N \yi, . . . ,y m ) is the m-point correlation function 

n^ N \ yi ,...,y m )= — — [ ■■■! V(y)dy m+1 ...dy N . (1.13) 

{N-m)\J R+ J R+ 

Our next main result shows the universality of the correlation kernel when M, N and 
Ni — > 00. As a polynomial in £, the algebraic equation (11.71) admits a unique solution 
£i(z) that is analytic in the upper half plane and vanishes at z — 00. If we let k 2 = 2, 
k 3 = 4 for a > 1 and k 2 = 4, /c 3 = 2 for a < 1, then at the points A^ .-i and X k , j = 2, 3, 
the root £i(z) will coincide with another root £j(z). Let 9j(z) be the following, 



9i( z ) = / €i( x )dx, 2 (z) = / ^ 2 (x)dx, 9 3 (z) = / ^ 3 (^)^, 

J\ 4 Jx k2 Jx k3 (1-14) 

k 2 = 2, k 3 = 4, a > 1, k 2 = 4, k 3 = 2, a < 1. 
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where the integration path is taken in the upper half plane. Then we can define a rescaled 
kernel Km,n(x, y) in a neighborhood of [Afc._i, A*-]. 



M-N 



K M>N (x,y) = (xy- 1 )^ e-( 9 ^^ + ^+^- 9 ^^- e -+^^ M ,7v(x,y), (1.15) 

where 0\,+ {x) and Oj t+ (x) are the boundary values of 9\(x) and 6j(x) on the positive side 
of the real axis. 

The rescaling from K M;N (x, y) to K MtN (x, y) in (11.151) does not affect the determinantal 
formula ( 11. 12ft . That is, we have 

K^ N \y u ...,y m ) = det (k MtN ( yj , y k j) (1.16) 

We then have the following universality result for the kernel Km,n{x,v)- 

Theorem 3. Suppose A in Theorem^ is positive. Let p(z) be the density function of F 
in Theorem^ Then for any x$ G (Aj., A2) U (A3, A4) and u, v G R ; we have 

,. 1 T y ( u v \ smir(u — v) 

hm — — — K K MN x + — — — -, x + -r-j—, — r = -, 7—. (1.17) 

N,M^ooMp(x ) \ Mp(x ) Mp(x )J ir(u — v) 

On the other hand, let x = Xk, k — 1, . . . , 4, then for u, v G R, we have 

1 T y ( u v 
Ar hm 2 K M,N x H r , x H a" 

' 00 (Afp^) 3 V (M Pfc )^ (M Pfe )V (L18) 

Ah>)Ai>) - Ai'MAifo) 

M — f 

where Ai(z) is the Airy function and p^, k = 1, . . . , 4 are the constants in 

Recall that the Airy function is the unique solution to the differential equation v" = zv 
that has the following asymptotic behavior as z — > 00 in the sector — 7r + e < arg(z) < n — e, 
for any e > 0. 

Ai(z) = — ] -^e~'^ (l + 0(z-*j) , -7r + e < arg(z) < vr - e, z -> 00. (1.19) 

3 

where the branch cut of 22 in the above is chosen to be the negative real axis. 

Since the limiting kernel Km,n{%-iD) takes the form of the Airy kernel (I1.18p . by a 
well-known result of Tracy and Widom [3U] , we have the following 



Theorem 4. Let y\ be the largest eigenvalue of B^, then we have 



Jim P ( ( Vl - A 4 ) (Mp 4 )^ <s) = TW{s), (1.20) 



where TW(s) is the Tracy-Widom distribution 

TW{s) = exp y- {t- s)q 2 {t)d?j , 

and q{s) is the solution of Painleve II equation 

q"(s) = sq(s) + 2q 3 (s), 

with the following asymptotic behavior as s — > oo. 

q(s) ~ — Ai(s), s — > +00. 

Remark 2. The results obtained in this paper are derived from the asymptotics of multiple 
Laguerre polynomials (See Section\^for details). In \2l^ , strong asymptotics of the multiple 
Laguerre polynomials was also obtained for the case when 1 — c = 0{M~ X ) and (3 — |. 
This corresponds to the case when M — N is finite. In [27^ , it was shown that, in that 
case, the origin will become an edge of the spectrum and the asymptotics of the multiple 
Laguerre polynomials are described by the Bessel functions near the origin. The statistical 
implication of |F7| / is that, when M — N is finite, the correlation kernel near the origin will 
be given by the Bessel kernel (See, e.g. [31]) instead of the Airy kernel U.18\) . This is an 
analogue to the situation when £ a? = I and M — N is finite, where the correlation kernel 
at the origin is also described by the Bessel kernel (See J73]/ and f3Ej). 

Until recently, most of the universality results for the Wishart distribution was obtained 
when the covariance matrix Ejv is the identity matrix [20], [21], [25] and [26]. More recent 
studies have extended these results to the spiked model proposed by Johnstone [2S], in 
which Ejv is a finite perturbation of the identity matrix [BJ, [7j, [E], P3]> [S3], [31] • This 
is the first few cases when universality results was obtained for a covariance matrix that 
is not a finite perturbation of the identity matrix. (See also [19] in which a different class 
of was studied) For theoretical reasons, the model studied in this paper is crucial in 
understanding of the phase transition behavior that occurs in Wishart ensembles. (See 
[7]). When the 2 intervals in the support of F closes up, a phase transition takes place and 
the correlation kernel will demonstrate a different behavior at the point where the support 
closes up. With the Riemann-Hilbert technique used in this paper, such behavior can be 
studied rigorously as in [12] (See also [27]). We plan to study this phenomenon in a further 
publication. For practical reasons, many covariance matrices appearing in fields of science 
are not finite perturbations of the identity matrix. In fact, covariance matrices that have 
groups of distinct eigenvalues are accepted models in various areas such as the correlation 
of genes in microarray analysis and the correlation of the returns of stocks in finance. 
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2 Multiple Laguerre polynomials and the correlation 
kernel 



The main tool in our analysis involves the use of multiple orthogonal polynomials and 
the Riemann- Hilbert problem associated with them. In this section we shall recall the 
results in [13] and [15] and express the correlation kernel K MjN (x,y) in (11.121) in terms of 
the multiple Laguerre polynomials. In Section HJ we will apply Riemann-Hilbert analysis 
to obtain the asymptotics of these multiple Laguerre polynomials and use them to prove 
Theorem [3j 

We shall not define the multiple Laguerre polynomials in the most general setting, but 
only define the ones that are relevant to our set up. Readers who are interested in the 
theory of multiple orthogonal polynomials can consult the papers [T], [2], [13], [22]- Let 
L ni ,n 2 ( x ) be the monic polynomial such that 



oo 

i+M-N -Mx 



L nin Jx)x l+M -"e- Mx dx = 0, % = 0, . . . ,m 



ni,n2 

oo 

T 

m,n2 



L ni n 2 (x)x i+M ~ N e~ Ma ' lx dx = 0, I = 0, ... ,712 - 1. 



(2.1) 



and let Q ni ,n 2 { x ) be a function of the form 

Qni,n 2 ( x ) — ^■n 1 ,n 2 ( X ) e ^ + ^n b n 2 ( X ) e ^ X > (2-2) 

where A l ni Tl2 (x) and A® in2 (x) are polynomials of degrees n\ — 1 and n 2 — I respectively, 
and that Q ni ,n 2 ( x ) satisfies the following 



oo 



x Qrii ,no (x\x dx \ i ■ i (2-3) 

' y 1, i = rii + n 2 — 1. 

The polynomial L„ i n2 (x) is called the multiple Laguerre polynomial of type II and the 
polynomials A\ n2 (x) and A a nx ri2 (x) are called multiple Laguerre polynomials of type I 
(with respect to the weights x M ~ N e~ Mx and X M ~ N e ~ Ma a; ) pQ^ [2]. These polynomials 
exist and are unique. Moreover, they admit integral expressions 
Let us define the constants hn},n 2 

and /irii,?i2 

to be 



ni,n 2 



/oo 
Ln-i,n 2 i x ') x & dx, 

,00 (2-4) 

*SU= / L nijn2 (x)x^ +M - N e- M ^dx. 
Jo 

Then the following result in [13] and [15] allows us to express the correlation kernel in 
(11.121) in terms of a finite sum of the multiple Laguerre polynomials. 
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Proposition 1. The correlation kernel in Km,n{x,u) $UB) can be expressed in terms of 
multiple Laguerre polynomials as follows 



i v — 1V1 

(xy) 2 (x -y)K MtN (x,y) = L NotNl (x)Q No>Nl (x) 



h {1) 

L No _ hNl {x)lJ No+hNl {x) 

n>N -i,Ni y Z - b > 



h (2) 

— ^N ,Ni-l \ x ) VJVcJVi+l [X) 



h {2) 



where N = N — N±. 

This result allows us to compute the limiting kernel once we obtain the asymptotics 
for the multiple Laguerre polynomials. In this paper, we will use the Riemann-Hilbert 
method to obtain such asymptotics and use them to compute the limiting kernel and 
to prove Theorem 2. The Riemann-Hilbert analysis used in this paper involves a 3 x 3 
Riemann-Hilbert problem and the analysis is similar to that in [TT] where random matrices 
with external source was studied (See also [27]). In what follows, we will restrict ourselves 
to the case when the limiting eigenvalue distribution is supported on 2 disjoint intervals. 
The case when the limiting distribution is supported on 1 interval will be considered in a 
separate publication. 



3 Stieltjes transform of the eigenvalue distribution 

In order to study the asymptotics of the correlation kernel, we would need to know the 
asymptotic eigenvalue distribution of the Wishart ensemble (11.31) . Let F N (x) be the empir- 
ical distribution function (e.d.f) of the eigenvalues of Bn ( 11.21) . The asymptotic properties 
of Fn(x) can be found by studying its Stieltjes transform. 

The Stieltjes transform of a probability distribution function (p.d.f) G(x) is defined by 

r°° 1 

m G {z) = / dG{x), z € C + = {z e C : Im(z) > 0} . (3.1) 

Given the Stieljes transform, the p.d.f can be found by the inversion formula 

1 f b 

G([a,b}) = - lim / Im {m G (s + ie)) ds. (3.2) 

e-,0+ J a 

The properties of the Stieltjes transform of Fn(x) has been studied by Bai, Silverstein and 
Choi in the papers [H], [3], pi], [5], |29j. We will now summarize the results that we need 
from these papers. 
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First let us denote the e.d.f of the eigenvalues of by if at, that is, we have 



dH N {x) = -^XX-- 



Furthermore, we assume that as X —>■ oo, the distribution weakly converges to a dis- 
tribution function H . Then as X — > oo, the e.d.f F^[x) converges weakly to a nonrandom 
p.d.f F, and that the Stieltjes transform m F of F(x) satisfies the following equation [T4"] . 

ESI 



= / t(l-c-lm F )-/ m - (3 - 3) 



Let us now consider the closely related matrix B_ N 

B N = ^XS^Xt. (3.4) 

The matrix B_ N has the same eigenvalues as together with an addition M — N zero 
eigenvalues. Therefore the e.d.f F N of the eigenvalues of B N are related to Fn by 

X 

F N = (1 - Cjv)-T[0,oo) + c N F N , c N = — . (3.5) 

where I[o,oo) is the step function that is on R_ and 1 on IR + . In particular, as X — > oo, 
the distribution F_ N converges weakly to a p.d.f F_ that is related to F by 

F = (1 - c)/ [0>oo) + cF (3.6) 

and their Stieltjes transforms are related by 

1 - c 

m_p(z) = hcmp^). (3.7) 

Then from (13.31) . we see that the Stieltjes transform mp(z) satisfies the following equation 



tdH(t) 



-i 



m E.( z ) = -\ z - c I Y+u^t ' ' (3 - 8) 



This equation has an inverse [I], [5] 



1 ft 

z (m F ) = + c / dff(t). (3.9) 

- m F J R 1 + tm F 
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3.1 Riemann surface and the Stieltjes transform 

We will now restrict ourselves to the case when the matrix has 2 distinct eigenvalues 
only. Without lost of generality, we will assume that one of these values is 1 and the other 
one is a > 0. Let < /3 < 1, we will assume that as N — > oo, Ni of the eigenvalues take the 
value a and N = N — Ni of the eigenvalues are 1 and that ^ — > (3. That is, as N — > oo, 
the e.d.f H^{x) converges to the following 

dH N (x) -> dH(x) = (1 - p)5i + /35 a . (3.10) 

By substituting this back into the (I3.9p . we see that the Stieltjes transform £(z) = mp_{z) 
is a solution to the following algebraic equation 

s = -! + c TT? + c TT<' (3 ' u) 

Rearranging the terms, we see that £ = itif_ solves the following 

za^ + (A 2 z + B 2 )i 2 + (z + BJt + 1 = 0, 
A 2 = (l + a), B 2 = a(l-c), (3.12) 
B x = 1 -c(l-/3) +a(l -c(3). 

This defines a Riemann surface £ as a 3-folded cover of the complex plane. 

By solving the cubic equation (I3.12p or by analyzing the asymptotic behavior of the 
equation as z — * oo, we see that these solutions have the following behavior as z — ■> oo. 

= -- + 0(z~ 2 ), 2-^00, 

£ 2 (^) = -l + c(1 ~^ +0(z- 2 ) } z^oo, (3.13) 

£ 3 (^) = -- + - + 0(z- 2 ), z^oo. 
a z 

On the other hand, as z — > 0, the 3 branches of £(z) behave as follows 



c 



=i?! + 0(z), Z^O, 
£ 7 (*) =-Ra + 0(*), z^O. 



(3.14) 



where the order of the indices a, f3 and 7 does not necessarily coincides with the ones in 
( 13.131) (i.e. we do not necessarily have a = 1, (3 = 2 and 7 = 3). The constants Ri and R 2 
are the two roots of the quadratic equation 

a(l - c)x 2 + (1 - c(l - /3) + a(l - c/3))x + 1 = 0. (3.15) 
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The discriminant D 2 of (13.151) is given by 

D 2 = (1 - c(l -/3) + a(l - c(3)f - 4o(l - c), 



(3.16) 



for a < 1, it is an strictly increasing function in (3 and hence 

D 2 > (1 - c + a) 2 - 4a(l - c) = (1 - c - a) 2 > 0, 

and for a > 1, it is a strictly decreasing function in (3 and hence 

D 2 > (1 + a(l - c)) 2 - 4a(l - c) = (1 - a(l - c)) 2 > 0. 

Therefore both R± and i?2 are real as they should be. 

The functions £j(z) will not be analytic at the branch points of £ and they will be 
discontinuous across the branch cuts joining these branch points. Moreover, from (13.141) . 
one of the functions £j(z) will have a simple pole at z — 0. Apart from these singularities, 
however, the functions £j(z) are analytic. 

3.2 Sheet structure of the Riemann surface 

In this section we will study the sheet structure of the Riemann surface C. As we shall 
see, the branch £1(2) turns out to be the Stieljes transform mp_{z) and its branch cut will 
become the support of F. 

3.2.1 The support of eigenvalues 

The branch points of the Riemann surface £ are the points on C in which ^ vanishes. 
They are also potential end points of supp(F) due to the following result in [T4"] . 

Lemma 1. [7^] If z £ supp(F) ; then m = mp^z) satisfies the following. 

1. m Gl\{0}; 

2. i supp(F); 

3. Let z be defined by $3. 9\) . then z'{m) > 0, where the prime denotes the derivative with 
respect to mp_ in $3.9\) . 

Conversely, if m satisfies 1-3, then z = z(m) supp(F). 

This lemma allows us to identify the complement of supp(F) by studying the real points 
m such that z'{m) > 0. 

Let us differentiate (13.111) to obtain an expression of ^ in terms of £. 

a 2 (l - c)£ 4 + 2(a 2 (l - c(3) + a(l - c(l - (3))^ 

(3.17) 

+ (1 - c(l -P) + a 2 (l - c(3) + Aa)i 2 + 2(1 + a)£ + 1 J . 
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dz 
di 



e(i + o 2 (i + <) 2 



In particular, the values of £ at ^ = correspond to the roots of the quartic equation 

a\l - c)£ 4 + 2(a 2 (l - c/3) + a(l - c(l - /?))£ 3 

+ (1 - c(l - /?) + a 2 (l - c/3) + 4a)£ 2 + 2(1 + a)£ + 1 = 1 ' ' 

Let A be the discriminant of this quartic polynomial, then when A > 0, the equation 
(13.181) will have 4 distinct real roots 71 < . . . < 74. Since < c, (3 < 1, the coefficients of 
(13.181) are all positive and hence all 7^ are negative. 
Let Afc be the corresponding points in the z-plane 

X k = -- + c^- + c-^—, fc=l,...,4. (3.19) 
Ik 1 + Ik 1 + aik 

Note that, from the behavior of z(£) in (13. lip , we see that near the points —1 and — , the 
function z(£) behaves as 

A 1 (3-20) 

1 + a£ a 

and hence i s negative near these points. From this and ( 13. 17ft . we see that z'(^) > on 
the intervals (—00, 7^, [72, 73], [74, 0) and (0, 00) and none of the points — 1 or — - belongs 
to these intervals. On [71,72] and [73,74], is negative. 

The images of these intervals under the map (13.111) then give us the complement of 
supp(F) in the z-plane. Let us study these images 

Lemma 2. The intervals (—00, 71], [72, 73], [74, 0) and (0, 00) are mapped by z(£) to (0, Ai], 
[A 2 , A 3 ], [A 4 , 00) and (— 00, 0) respectively. Furthermore, we have A 2 < A3. 

Proof. Since none of the points — 1, — - and belongs to these intervals both z(£) and 
are continuous on these intervals. Moreover, z(£) is strictly increasing on these intervals. 
Therefore the images of these intervals are given by 

2 ((-00,71]) = (z(-oo),z(jx)] = (0, Ai] 
z ([72, 7s]) = [2(72), 2(73)] = [A 2 ,A 3 ] 
2 ([74,0)) = [z(74),*(0-)) = [A 4 ,oo) 
2 ((0,oo)) = (*(0+),*(oo)) = (-oo,0). 

where the ± superscripts in the above indicates that the function is evaluated at ±e for 
e — > 0. Finally, since 73 > 72, we see that A3 > A2. This concludes the proof. □ 

Therefore the complement of supp(F) is given by (recall that F has a point mass at 0) 
supp(F) c = (-00, 0) U (0, Ai) U (A 2 , A 3 ) U (A 4 , 00). (3.21) 
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Thus if Ai < A2 and A3 < A4, the support of F will consist of the 2 intervals [Ai, A2] and 
[A3, A4]. We would like to show that whenever A > 0, we have Ai < A2 < A3 < A4. To do 
this, let us take a look at the zeros of the function from the point of view of branch 
points. 

The Afc are the ^-coordinates of the zeros of on C We can treat (I3.12p as a 
polynomial in £ then A& will be the zeroes of its discriminant D 3 (z) = (az) IIt<j(£» — £j) 2 - 

D 3 (z) = (1 - a) V + {2A 2 2 B 1 + 2A 2 B 2 - 4A 3 2 - \2aB x + 18aA 2 )z 3 

+ {B\ + A\B\ + 4A 2 B 1 B 2 - Y1A\B 2 - YlaB\ + 18a5 2 + \%aA 2 B x - 27a 2 )z 2 
+ {2B X B 2 2 + 2A 2 B 2 B\ - \2A 2 B 2 2 - 4Bfa + IZaB^z + B\B 2 2 - 4Bf . 

(3.22) 

The zeros of (13.221) also correspond to the branch points of the Riemann surface C. These 
branch points are given on C by (A&, 7*.), for k — 1, . . . , 4. 
Let us rename the Afc as z^. 

{Ai, . . . , A 4 } = {zi, . . . , Z4} , 

where the above equality is only an equality between the sets. The ordering of the points 
does not necessarily coincide. We shall order the points Zj such that z\ < z 2 < z 3 < z 4 . 

Let J be the union of the intervals [^1,^2] an d [-23,-24]. Since the leading coefficient of 
D 3 (z) is (1 — a) 2 > 0, we see that the sign of D 3 (z) and hence the 3 roots of the cubic 
polynomial (13.121) behave as follows for zEl. 

1. z Gl\ J, D 3 (z) > =>■ f has 3 distinct real roots 

(3.23) 

2. z G J, D 3 (z) < =^> £ has 1 real and 2 complex roots. 

In particular, since the roots coincide at the branch points, the jj are the values of the 
double roots of the cubic (I3.12p when z = \j. We then have the following lemma. 

Lemma 3. If a > 1, then —1 G [71,72] and — i G [73,74]- On the other hand, if a < 1, 
then —- G [71,72] and —1 G [73,74]- This means that z(£) has a singularity on both of the 
intervals [71,72] and [73,74]. 

Proof. We will show that cannot be continuous on either [71, 72] or [73, 74]. This means 
that both intervals must contain a singularity of z(£), which is one of the points —1, — . 
This would imply the lemma. 

Let us assume that none of the points —1 or —1 belongs to [71,72]- Then the function 
z(£) is continuous on [71,72]. Moreover, from the behavior of z at these points (13.201) . 
we see that none of these points belongs to [72,73] either and hence z(£) is continuous on 
[71 j 73]- We also seen in the proof of Lemma [2] that A2 < A3. Let us also note that Ai < A4. 
If not, then (13.211) would imply that the support of F_ is empty, which cannot happen. 
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Therefore there are only 3 possible ways of ordering the points Ai, A2 and A3. 

1. A 3 > Ai > A 2 ; 

2. Ax > A 3 > A 2 ; (3.24) 

3. A 3 > A 2 > Ai. 

We will show that z(£) cannot be continuous on [71,73] in any of these cases. 

First note that by the remark after (13.201) . we see that is negative between 7! and 
72 and hence we must have Ai > A2 if z(£) is continuous on [71,72]. We therefore rule out 
the third case in (13.241) . 

Let us now assume A3 > Ai > A2. Then by (13.231) and the fact that Ai < A4, we see 
that for z G [A2, Ai], there is only one real m such that z = z{m). Let zq G [A2, Ai]. Since 
z{£) is continuous between 71 and 72, the interval [A2, Ai] is in the image of [71,72] under 
z(£). Hence there is at least one point m in [71,72] such that z(m ) = z . Similarly, since 
z(£) is continuous between 72 and 73, the interval [A 2 , A3] is in the image of [72,73] under 
z(£). Since A3 > Ai > A2, the interval [Ai,A2] is a subset of [A2,Aa] and hence there is at 
least one point mi G [72,73] such that z(m\) = z . This contradicts the fact that there is 
only one real m such that zq = z(m) in [A2, Ai]. 

Let us now consider the case when Ai > A3 > A2. In this case, for zq G [A2, A3], there is 
only one real m such that z = z(m). Again, by continuity of z(£) on [72,73] and [71,72], 
we see that there is at least one point m in [72,73] and another point mi G [71,72] such 
that z(mi) = z(mo) = Zq. This leads to a contradiction and hence this cannot be the case 
either. 

By using the same argument, we can show that z(£) cannot be continuous on [72,74] 
and hence one of the points —1, — ^ must be in [71,72] while the other one is in [73,74]. 
This implies the lemma. □ 

We can now show that the points A& are ordered as Ai < A2 < A3 < A4. 

Lemma 4. The A& satisfies Ai < A2 < A 3 < A4. 

Proof. As we have seen in the proof of Lemma El there are only 3 possible ways of ordering 
the points Ai, A2 and A3, which is indicated in (13.241) . We will show that the cases 1 and 
2 in (13.241) are not possible. 

Let us assume A3 > Ai > A2. By LemmaEJ there is a singularity of z(£) in [71,72]. Let 
us call this singularity s . Let z G [A 2 , Ai]. Then there is only one real m G ffi\{0, —1, — ^} 
such that z(m) = z . By the continuity of z(£) on [72, 73], there exists a point m G [72, 73] 
such that z(m ) = z . On the other hand, let R > Ai, since z(£) is continuous on (s,7 2 ] 
and z(s + e) — > +00 for small e > 0, there exists 6 > such that [A2, -R] is contained in 
the image of [so + 5, 72] under the map z(£). In particular, there exists another mi on 
[sq + 5,72] such that z{m\) = z . This leads to a contradiction. 

Let us now assume Ai > A3 > A 2 . Then for zo G [A 2 ,A3], there can only be one real 
m G K\ {0, —1, — -} such that z(m) = Zq. Applying the continuity argument, we again see 
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that there is an mo G [72,73] such that 2 (mo) = z and that for any R > A3, there exists 
5 > such that [A2, -R] is contained in the image of [sq + 5, 72] and hence there is another 
point mi G [sq + 5, 72] such that z(mx) = Zq. This again leads to a contradiction and hence 
we must have A 3 > A 2 > Ai. 

By carrying out the same argument for the points A2, A3 and A4, we see that the only 
possible ordering of these points is A 2 < A 3 < A 4 . Hence we must have Ai < A 2 < A 3 < 
A 4 . □ 

We therefore have the following condition for the support to consists of 2 intervals. 

Theorem 5. Let A be the discriminant of the quartic polynomial 

a 2 (l - c)£ 4 + 2(a 2 (l - c(3) + a(l - c(l - (3))^ , . 

+ (1 - c(l -/3) + a 2 (l - cf3) + 4a)£ 2 + 2(1 + a)£ + 1 = 0. 1 ' ' 

then the support of F_ consists of 2 disjoint intervals if and only if A > . 

Proof. By Lemma H] and f!3.2ip . we see that if A > 0, then the support of F_ consists of 
2 disjoint intervals. If A < 0, then there can be at most 3 distinct roots to the equation 
( 13.171) and hence there can be at most 3 end points to the support of F_. This means that 
F_ cannot be supported on 2 disjoint intervals. □ 

From this theorem, we see that all the A& are positive. 

Corollary 1. The A/, satisfies < Ai < . . . < A 4 . 

Proof. From Theorem El we see that the limiting eigenvalue distribution is supported on 
[Ai,A 2 ] U [A3,A 4 ]. Therefore all A& > 0. We want to show that Ai > 0. Since the Afc are 
solutions to (13.2211 . Ai = if and only if the constant term of ( 13.221) is zero. This constant 
term is 

B 2 2 (Bj - AB 2 ) = a 2 (l - cfD 2 > 
where _D 2 is given by ( 13 . 1 6D . This concludes the proof. □ 

Let us now compute the probability density F_ when it is supported on 2 disjoint inter- 
vals. 

Theorem 6. Suppose A > 0. Then the p.d.f F_ is supported on [Ai, A 2 ] U [A3, A 4 ] with the 
following density dF(z) = p(z)dz 



3 

2^ 



(3.26) 
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where D%(z) and r(z) are given by 



4 

D 3 (z) = (l-a) 2 H(z-\ j ), 
r(2) = — I *-z + 5 U + —5- H 5 »— I z 




-1 




T/ie cu&e rooi in $3.26}) is chosen such that \YA G R /or i G R and the square root is 
chosen such that y/A > for A > 0. 

Proof. Let z £ [Ai, A 2 ] U [A3, A 4 ] be a point in the support on F. Then as a polynomial in 
£, (I3.12p has the following solutions at z = z . 

£ R (z ) = -l(** + — J + ^7 + ^r, 

3 \ a a2; J 

I (A 2 B 2 \ _ , 2 



where w± are 



6^0) = " o — + — + ^v^T- + w ( 3 -27) 
3 \ a az / 

62(^0) = — 5 ( — + — ) +^ 2 ^M7 + ^^r~, 
3 V a az J 



u± = W^ftj (3 ,S) 



2iri 



and a; = e~3~ is the cube root of unity. The functions r(z), 1)3(2) and constants v4 2 , i? 2 
defined in the statement of the theorem. The branches of the cube root in (13.271) and the 
square root in (13.281) can be chosen arbitrarily, as long as all the branches are the same. We 
will not treat ( 13.27P as analytic functions in zq, but merely consider them as the different 
values of the roots of ( 13.121) at the point z = z . 

Since D 3 (z) < on supp(F), u±(z) G R on the support of F. Let us choose the square 
root and cubic root in (13.271) in the way that is indicated in the statement of the theorem. 
Then for any point z G supp(F), (,r(z ) is real while Ch( z o) and (,i 2 (z ) are complex 
conjugate to each other. Since the analytic continuation of rriF_{z) on the positive side of R 
becomes one of the complex roots on supp(F), we see that the imaginary part of m^(zo) 
must coincide with the imaginary part of either ^(^o) or £,h( z o)- I n particular, we have 

|Im(m z (*o))| = \ l&i(*o) -&.(*>)! • (3.29) 
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From the inversion formula and the fact that the p.d.f. F_ is non-negative, we see that the 
probability density for F_ is given by 

dF = ^- \£ h {z) - £ h (z)\dz, z e supp(F). 

Z7T 

By substituting ( 13.271) back into this equation, we arrived at the theorem. □ 



3.2.2 Branch cuts of Riemann surface 

In this section we will show that the solution £,i(z) (13.131) of (13. 12ft coincides with the 
Stieltjes transform mp_(z) in C + U IBL Furthermore, when the support of F_ consists of 2 
disjoint intervals, the function £1(2) and hence mp_(z) will not be analytic at any of the 
points Afe, k = 1, . . . , 4. 

The solutions £j(z) of (13.121) will not be analytic at the branch point (Xkj^k) if £j(Afc) = 
7fc. Since all the A^ are on the real axis, and that the only possible pole of these functions 
is at z = 0, there exists analytic continuations of the £j(z) in C + that are continuous up 
to K. \ {0}. We have the following lemma 

Lemma 5. Let £j(z) be the solutions of $3.12\) in (C + Ul)\ {0} that has the asymptotic 
behavior $3.13\) . Then for j = 1,2,3, £j(Ai) = 71 if and only if £j(\2) = 72- Similarly, 
&C*s) = 73 if and only i/^(A 4 ) = 74- 

Proof. Suppose £j(Ai) = 71 for some j. Then takes the value of the double root of (I3.12p 
at Ai and by continuity, £j becomes one of the complex roots of (13.121) on the interval 
[Ai, A2]. This implies that at A2, £7 will take the value of the double root of (13. 121) again. 
Therefore we have £j(A2) = 72- Similarly, we see that if Cj(^2) = 72, then we also have 
<v/(Ai) = 7i- Hence £j(Ai) = 71 if and only if ^j(A 2 ) = 72- Applying the same argument to 
the points A 3 and A 4 , the lemma is proven. □ 

Let us now show that the function £i(z) is in fact the Stieltjes transform mp_(z). From 
the asymptotic behavior of the £j(z) (13. 131) and the fact that the Stieltjes transform m^(z) 
solves (13. 121) and vanishes as z — > 00, we see that 

mpiz) = ^(z), zE (C + UM) \ {0}. (3.30) 

For c < 1, it was shown in [14J that F has a continuous density and hence the Stieltjes 
transform mp(z) does not have any poles. Therefore, by ( \3.7L we see that m^(z), and 
hence £1(2), has the following singularity at z — 0. 

£i(*) = -i^ + 0(l), z^0. (3.31) 
We will now show that for all A^, we have £i(Afc) — lu- 
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Proposition 2. Let £1(2) be the solution of ^3.12^ in (C + Ul)\ {0} with the asymptotic 
behavior indicated as in $3.13\) . Then we have £i(Afc) = 7fc for k — 1, ... ,4. In particular, 
£1(2) is not analytic at any of the points A&. 

Proof. Suppose for some A& we have £i(Afc) 7^ 7^ and that the support of F_ is on the left 
hand side of A&. Then we have z' (£i(Afc)) 7^ 0. By lemma[U we must have z' (mp^Xk + 5)) > 
for small enough 5 > as A& + 5 does not belong to the support of F and that rrip is real 
on Afc + 5. Then by continuity, we have z' [mpiX^j) > 0. (As we assume z' [rapiX^j) 0) 
Since £1(2) does not coincide with the double root of (I3.12p at A&, it must be real in 
[Afc — e, Afc + e] for some small e > 0. For small enough e, we can assume that neither 0, -1 
or — - are inside this interval. 

a 

Therefore in the interval [A& — e, A& + e], the Stieltjes transform m,F_(z) is real and by 
continuity, z' (m^) > inside this interval. Since this interval consists of points that 
belong to the support of F_, this would imply that for some point z G supp(F), we have 
z' (mp) > and mp 6R\{0, —1, — ^}. This contradicts lemma [T] and hence £i(Afc) = 7&. 
By using exactly the same argument, we can prove the proposition for A& when the support 
lies on the right hand side of A&. □ 

This proposition implies that £1(2) is not analytic at the points A^ for k = 1, ... ,4. 
With this information, we can now find out the structure of the Riemann surface L. 

Proposition 3. If a > 1, then £3(2) has branch points at A3 and A4 and ^{z) has branch 
points at X\ and A2. On the other hand, if a < I, then £3(2) has branch points at Ai and 
X2 while £2(2) has branch points at A3 and A4. 

Proof. We will prove the statement for the case a > 1, the case a < 1 can be proven by 
similar argument. Suppose a > 1, then for large enough zo > A4, all the functions £j are 
real for z G M and we have £1(^0) > £3(^0) > £2(^0) by (13.131) . This ordering must preserve 
at A4 as the roots cannot coincide between zq and A4. At A4, one of the roots must coincide 
with £x and from the ordering £i(zo) > £a( 2: o) > £2(^0), we must have £i(A 4 ) = £s(A 4 ) = 74. 
By lemmaEl we also have £3^3) = 73. On the other hand, for small enough Ai > e > 0, all 
three £j will be real. From the asymptotic behavior of £1(2) at (13.311) we see that for small 
enough e, £i(e) will be smaller than both £2(6) and £3(6). From the asymptotic behavior 
of £2 and £3 at z = —00 (13. 13[) and the fact that these 2 functions has no singularity and 
cannot coincide on (—00, Ai), we see that at z — e, we must have £3(6) > £2(6)- Therefore 
we have £3(6) > £2(6) > £i(e)- This ordering must again be preserved at Ai. Therefore we 
must have £i(Ai) = £2(Ai) = 71. By lemma El we also have £2^2) = 72- This shows that 
for a > 1, the branch points of £3 are A3 and A4 while the branch points of £2 are Ai and 
A2. The proof for the case a < 1 follows from similar argument. □ 

We will now define the branch cuts of the function £i(z) to be [Ai,A2] U [A3, A4] and 
the branch cut for £2(2), £3(2) to be [Ai, A2] and [A3, A4] respectively for a > 1 and [A3, A4], 
[Ai, A2] respectively for a < 1. We then have the following relations between the £j on the 
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Figure 1: The branch cut structure of the Riemann surface C when a > 1. 



branch cuts 



1. a > 1, 

2. a<l, 



&, T (z), z e [Ai, A 2 ]; 
£3,=f0), 2 e [A 3 , Aj. 

6,=f(*)> 2 e t A 3, A 4 ]; 



(3.32) 



where £^±(2) indicates the boundary values of £j on the ± sides of the branch cuts. 
The branch cut structure of the Riemann surface £ is indicated in Figure [H 
We will now define the functions 9j(z) to be the the integrals of £j(z). 



e 1 (z) 



£i(x)dx, 9 2 {z) 



A 4 



£, 2 (x)dx, 9 3 {z) 



^(x)dx, 



(3.33) 



a > 1, k 2 



a < 1. 



where the integration paths of the above integrals are chosen such that they do not intersect 
the real axis, except perhaps at the end points. 

Then from (13.131) . (13. 14|) and (13.311) . we see that the integrals ( I3.33P have the following 
behavior at z = 00 and z = 0. 



9i(z 
9i(z 
9 2 {z 
9 2 {z 

9 3 (z 
9 3 (z 



00, 



-logz + h + (z- 1 ) , 1 
-(1 -c) log* + 0(1), z^O, 
-z + c(l-0) log z + l 2 + 0(z- 1 ) 
0(1), z^O, 

+ c(3 log z + l 3 + O (z' 1 ) , z 



00, 



(3.34) 



z 
a 

0(1) 



00, 



for some constants h, l 2 and ^3. 

From the behavior of £j(z) on the cuts, we have the following analyticity properties of 



the 9j(z). 
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Lemma 6. The integral 9i(z) is analytic on C\(— oo, A4] and continuous up to R\{0}. The 
integrals #2(2) and 6^{z) are analytic on C \ (—00, X k2 ] and C \ (—00, Afc 3 ] respectively and 
are continuous up to IR. Across the real axis, they have the following jump discontinuities. 



0i,± 




n 

= U j,T 


[z) 


+ 0i,±(A fcj .J, z G [X kj -i, A fcj J, j = 1, 2, 




» 


= K- 


[z] 


— 2cf37ii, z G [A 2 , A 3 ], a > 1, 






= K- 


[z] 


- 2c(l - /3)7ri, 2G[A 2 ,A 3 ], a < 1, 






= 0i,- 


[z] 


— 2C7TZ, Z G (0, Al], 






= K- 


(z) 


— 27ri, 2 G (—00, 0], 


02,+ 




= e 2t - 


[z] 


+ 2c[l-p)ni, z G (—00, Afc 2 _i], 


o 3 ,+ 






(2) 


+ 2cf3iri, z G (-00, Afc 3 _i]. 



(3.35) 



Proof. From the jump discontinuities of £7(2) in (13.321) . we see that 



^i(x)c?a; 



0j#{z), z G [A fcj _i, A fcj ], j = 2,3, 



By comparing this and (13.331) . we obtain the first equation in (13.35!) . 
Let us now show that 

02,+ (*) = 02,-(*) + 2c(l - /3)7Ti, Z G (-OO, Afca-i]. 



(3.36) 



The corresponding equation for 9\{z) and 9^{z) can be proven in a similar way. 

Let z G (—00, Afe 2 _i] and let T± be contours from X k2 to 2; such that T± are in the upper 
and lower half planes respectively. Then we have 



02,±(z) = / £2(3:) dx. 
Let T be the close loop on C such that T = T + — T_, then we have 

02,+{z) = 02,-(z) + j> &{x)dx, z G (-00, Afc 2 _l]. 
By Cauchy's theorem, we can deform the loop T such that T becomes a loop around z = 00. 



By computing the residue, we arrived at (13.36!) . 

On the other hand, we can also deform T so that it becomes a loop enclosing the branch 
cut [Afca-i, Afc 2 ] of ^(z). Then we have 



^ 2 (x)dx 



H: 2 



(3.37) 



From (13.321) . we have ^2,-(z) = £1,4.(2;). Since £ 2 (z) and £1(2) are the 2 complex roots of 
(13.121) on [Afc 2 _i, Afc 2 ], we have 

£,2,+( x ) ~ £2 ,-( x ) = fi-fa) ~ £iA x ) = -2ilm (fi,+(a:)) , x G [A fc2 _i, A fc2 ]. 
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By comparing ( 13.371) and ( 13.361) . we obtain 

lm(^+(x))dx = c(l -P)n. (3.38) 



By using similar argument, we see that 

03,+(z) = 03,-(z) + 2c{3iri, z G (-oo, A fc3 _i] 



/ Im (£i +(x)) dx — c/37r, 



(3.39) 



1>+ (*) = e x ,_(*) - 2cttz, ze(0,Ai], 
1)+ O) = 0i -0) - 2tt?, z G (-oo, 0]. 



From these, we have 



^k 9 r^k 3 



and hence 



Im (6,+ (a?)) = c(l - /3)vr, / Im (£i,+0)) = c/3vr. (3.40) 

fea-i 



/■A4 

1)+ (z) = 0i,-O) -2i / Im(^ li+ (x))da; 



A 3 

2cj37ii, a > 1; 

2c(l - /3)7rz, a < 1. 

for 2 G [A2 , A3] . This proves the lemma. □ 

From the proof of this lemma, we see that the mass of F is distributed between [Ai, A 2 ] 
and [A 3 , A 4 ] in the following way. 

Corollary 2. Let k?, = 2, k% = 4 for a > 1 and &2 = 4, k% = 2 for a < 1. Then the mass 
of the measure F_ is distributed between [Ai, A2] and [A3, A4] in the following way. 

F([A fca _i,A fca ]) = c(l-/3), F([A fc3 _ 1 ,A fc3 ])=c/?. (3.41) 

We will conclude this section with the following results on the relative sizes of the 
Re (9j(z)), which are essential in the implementation of the Riemann-Hilbert method. 

Lemma 7. The real parts of 9j{z) are continuous across R. Let us denote the real parts 
by Ke(8j(z)). Then we have the following 

Re(0i(*)) - Re(0i(Afc,.)) > Re(0;(s)), x eR + \ [X kj - U \ kj ], 3 = 1, 2. (3.42) 
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Proof. First note that, by ( 13.321) . we see that £j,±(z) are complex conjugations on [Ai, A2] U 
[A3, A4]. Therefore the real part Re(£j(z)) of £j(z) and hence Re(0j(z)) is continuous across 
the real axis. 

From the proof of Proposition [31 we have, for j = 2, 3, 

£i0)>&0), ^g[a 4 ,oo), 

Since and £j(z) are complex conjugate solutions on [Afe._i,At.] and that £1(2) is the 
branch with positive imaginary part (recall that £i(z) = itlf_{z) and from (13.21) . we see that 
the imaginary part of mp_(z) is positive on supp(F)), we have, by (13.271) . the following 
behavior near the branch points 



ix{z) = jkj-i + iv{Jz- Xkj-i + O ((z - Xkj-i)) , z Afe^-i, 



ij{z) = jkj-i - iv{Jz- A fej _i + 0((z- A fcj _i)) , z -> Xkj-i, 



€i(z) = Tfcj + V 2\J Z - ^ + O ((z - X k .)) , z -> X kj , 

£,A z ) = ik 3 - v{ Jz - x kj +o((z- x kj )) , z -> x kj . 



(3.44) 



where v\ and are real constants and the branch cut of the square root in the first 2 equa- 
tions is chosen to be the positive real axis, with the branch chosen such that v\_\J z — X kj -\ 
is positive on the left hand side of the positive real axis; while the branch cut of the square 
root in the last 2 equations is chosen to be the negative real axis, with the branch chosen 
such that v 3 2 \fz ~ ^k 3 -i is positive on the positive real axis. From (I3.44p . we see that there 
exists 5 > such that 

Re (&(z)) > Re (0(z)) , z G (A fc „ A fe , + 6), 

Re(£i(z)) <Re(^(z)), z G (A fcj _x - 5, A fcj _i). l ' ' 

The ordering (13.451) must be preserved until z hits another branch cut. Therefore from 
(jS35D and (E5D , we have 

Re (6(z)) > Re , « G [A fe ., 00) \ supp(F), 

Re (&(z)) < Re (£,-(z)) , z 6 (0, A fci _i] \ supp(F). 

We will now show that these inequalities hold for the whole intervals (0, X kj -i\ and [Afe,., 00). 
Let I 7^ j. Suppose either of the following happens in [Afc,_i, Afe,], 

Re (&(z) - > 0, if X kl > X kj , 

Re(^(z)-^(z))<0, if X kl < X kj ■ [ ' ' 

Then by (13.461) . there must be a point z\ G [Afe,_i, Afe,] such that 

Re(0(^i)) = Re(£i(z!)), 

(3-48) 

dz 



-(Re(^i)-^i)))>0. 
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The values of the £(z) at the point Z\ are given by ( 13.271) with £j(z) = £r(z) while £1(2) is 
one of the complex roots Oil- 2 ) or £hi z )- Taking the difference between the real parts of 
£r(z) and ^(z) (or £i 2 {z) which has the same real part as Ch( z ))i we have 

Re (&(*)- £i(z)) =2(^17+^7-) (3.49) 



Since z\ is a point on supp(F), from Lemma [U we see that the derivative C'r( z i) °f the real 
root at zi must be non-positive. By ( 13.271) . this derivative is given by 



?r( z o) = o-4 + -r (^7+^1) . (3.50) 



1 B 2 d 

3 az\ dz 

By ( tHUHD and fl3~5T?|) . we see that if £ (Re (^(z x ) - £i(*i))) > 0, then (13301 will be positive, 
this contradicts Lemma [Hand hence Re (£j(z) — £1(2)) cannot vanish in [A^-i, A^J. This, 
together with ( 13.461) implies 



Re(£i(z)) > Re{Zj(z)) , ze[X kj ,oo), 
Re (6(2)) < Re ze (0,A fc ._J. 



(3.51; 



Therefore we have 



/ Re (£i(x) - 0( x )) > 0, z > A fej , 

Re (£i(a;) - £j(x)) dx > 0, < 2; < A fcj _i, 

Since £1(2) — £7(2) is purely imaginary on [Aj^._i, A&J, we have 

Re(£i(x)-^-(x))da; = 0, (3.52) 



kj -i 



Then from the definition (13.331) of the 0j(z), we obtain (13.421) . □ 

The final result in this section concerns about the behavior of these real parts in a 
neighborhood of the branch cuts. 

Lemma 8. The open interval (X^-i, Xk 3 ) , J ' = 2, 3 has a neighborhood Uj in the complex 
plane such that 

Re(9j(z)) + Re {9 1 (X k .)) > Re(6 1 (z)) > Re(0 { (*)) + Re (0i(A fc ,)) , (3.53) 
where I = 2, 3 and I 7^ j . 
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Proof. Since £1(2) — £3(2) * s purely imaginary on [Afc._i, AfcJ, we have 

/ Re fa(x) - ^(x)) rfx = Re (9 1 (z) - 9 x {\ k .) - 9j(z)) = 0, z G [A fc ,_ 1; A fc .], (3.54) 

On the positive and negative sides of (Afe._i, AfcJ, the derivative of the functions Q\,±{z) — 
0j,±(z) are given by £i,±(z) —£j,±(z) and are purely imaginary. In fact, since £i,+ (z) = mp, 
we see that £1,4.(2) ~ £j-( z ) = ^ip(z) where p(z) > is the density function of F. 
On the other hand, by the jump discontinuities (13.321) . we see that £i,_(z) — £j,+{z) = 
—2nip(z). Hence by the Cauchy Riemann equation, the real part of 9\{z) — 9j(z) — #i(Afc.) 
is decreasing as we move from (A&._i, A&-) into the upper half plane. From (13.541) . we see 
that Re(9i(z) — 9j(z)—9i(\k J )) < for z in the upper half plane near (A^-i, A^). Similarly, 
we also have Re(9 1 (z) — 9j(z) — #i(Afc.)) < for z in the lower half plane near (At.-i, At.). 
Therefore in a neighborhood of (A^ _i, Afc.), we have 

Re(9 1 (z)) < Re(9 j (z) + 0i(A fe ,)). (3.55) 

Now by Lemma d we see that, if I = 2, 3 and I 7^ j, then we have 

Re(9 1 {z)) > Re(9 l (z) + 9 1 (X h )), z G (A fe ,_ x , X kj ). (3.56) 

From (13.561) and (I3.55p . we see that (I3.53j) is true. □ 



4 Riemann-Hilbert analysis 

We can now implement the Riemann-Hilbert method to obtain the strong asymptotics for 
the multiple Laguerre polynomials introduced in Section [5] and use it to prove Theorem [3J 
The analysis is very similar to those in [TTJ (See also \27\). 

Let C(f) be the Cauchy transform of the function f(z) G L 2 (M + ) in M + 

C(f)(z) = I ^-ds, (4.1) 



2vri 

and let Wi(z) and w 2 (z) be the weights of the multiple Laguerre polynomials. 

Wl ( z ) = z M - N e~ Mz , w 2 (z) = z^-V^, (4.2) 
Denote by K\ and k 2 the constants 

Ki = -2m f/i£Li at,) , t 2 = -2m ( hf ] 



v N -l,N 1 J ' n 2 — y"N ,Ni-l 

Then due to the orthogonality condition (12. ip . the following matrix 

/ Pnq^Az) C{P NqjNi w 1 ){z) C(P NotNl w 2 )(z) 
Y(z) = KtP No -i lNl (z) KlC(P No -i jNl Wi)(z) KiC(P No -i jNl w 2 )(z) I (4.3) 
\«2-PjVo,JVi-l(«) /«2C(PAr ,7V 1 _iWi)(z) K 2 C{P N0)Nl _ X W 2 ){z) ^ 
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is the unique solution of the following Riemann-Hilbert problem. 

1. y{z) is analytic in C \ R + , 

'1 wi{z) w 2 {zf 

2. V+U) = V'_(:) | 1 | . ; e 

1 

z N 

3. Y{z) = (I + 0(z- 1 )) [ z- No 






(4.4) 



oo, 



4. y(z) = o(i), z-^o. 

On the other hand, the following matrix 

/ -27riC(A 1 No!Ni w 1 + A a NQ>Ni w 2 ) 2m 2mA a No<Ni \ 

X(z) = — h NNl C [Aj^ 0+ltNl Wi + A^ o+1Ni w 2 ) h N Nl Ax 0+ltNl h No ^ Ni A^ Q+1 ^ Ni 

\ ^N ,Ni^ (^A r ,Afi+l u 'l + ^Wo.iVi+l^) ^JV ,JVi^iVo,JVi+l ^0,^1^0,^1+1/ 

is the unique solution of the following Riemann-Hilbert problem 

1. X(z) is analytic in C \ M + , 

1 N 

2. X+(;z) = X_(s) I -w 1 (z) 1 I , 2 e 



f z- N 
3. X{z) = (l + 0(z~ 1 )) ( 

z Nl 



(4.5) 



4. X(z) = 0(l), z^O. 

In particular, by the uniqueness of the solution of the Riemann-Hilbert problems (14.51) and 
flOD , we see that 

X(z) = Y- T {z). (4.6) 

The proof of these statement can be found in [22] ■ By a similar computation as that in [TU] 
and |13| . we can express the kernel ( 12. 5p in terms of the solution of the Riemann-Hilbert 
problem Y(z). 

(xy)^ (e- M y [y+ l (y)y + (x)] 21 + e-^- 1 » [y+^y+^y 

= 2^y) ' 

(4-7) 

(xy) 5 

■ y) 

where Am and A™ are the 21 th and 31* h entries of A. 



2m\x — y) 



e 



My e -Afa ^ F- 1 (y)F + (a;) 
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4.1 First transformation of the Riemann-Hilbert problem 

We should now use the functions 9j(z) constructed in Section [3221 to deform the Riemann- 
Hilbert problem (14.41) . Our goal is to deform the Riemann-Hilbert problem so that it can be 
approximated by a Riemann-Hilbert problem that is explicitly solvable. Before we deform 
the Riemann-Hilbert problem, let us make the following observation. 

As pointed out in [8], Lemma [U applies to any distribution G whose Stieltjes transform 
itlg{z) (13.11) satisfies an equation of the form 



1 ft 

z{m G ) = hc G / — — dH G (t). 



for some constant c G and distribution H G {t). Hence if we replace the constants c by 
cn = jf, ft by Pn = an d H{t) by H N (t) = (1 — (3n)5i + Pn$o,i the results we obtained in 
Section [3] remain valid. In particular, since the discriminant A in Theorem [6] is continuous 
in the parameters c and P, if A > 0, there exists large enough M, N and N\ such that the 
discriminant A^r evaluated with the parameters and is positive. Let Bf and B 2 
the following constants 

B* = a(l - c N ), Bf = 1 - c N (l - p N ) + a(l - c N p N ). 

Then for large enough M, N and N ly the equation 

z< 3 + (A 2 z + B?)£ 2 + (z + B?)£ + 1 = 0, (4.8) 

will define a Riemann surface with 4 distinct real branch points < A^ < . . . < A4 . In 
particular, let the roots of this equation be ^(z) and define Of (z) as in (I3.33P 



0?{z)= tf(x)dx, 0»(z)= e 2 \x)dx, 0»(z)= $(x)dx, 

J *¥ ' J ^ 2 J ^ 3 (4.9) 

k 2 = 2, h = A, a > 1, k 2 = 4, h = 2, a < 1. 

Then all the results we obtained in Section f3. 2. 21 will apply to the 0f(z) and £j(z) with c 
and P replaced by cat and Pn- In particular, Qf(z) has the behavior near z — > 00 and z = 
as indicated in (13.341) with c, P replaced by cat and Pn and lj replaced by some constants 
l N 

We will now start deforming the Riemann-Hilbert problem (14.4j) . First let us define the 
functions gf{z) to be 

g»{z) = 0?(z) + (1 - c N ) logz, <^(z) = 0f (z) + 0f_(A£) + z, 
^)=^) + <-(A£) + ~. (410) 

CI 

where the branch cut of log 2 in gf{z) is chosen to be the negative real axis. 
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We then define T(z) to be 

T{z) = diag (e-^.e-^.e-^) Y(z)diag f e W« e^M.eWM 

The matrix T(z) will satisfy the following Riemann-Hilbert problem. 

1. T(z) is analytic in C \ R, 

2. T+(z) =T_(z)J T (z), zeR, 

3. T(z) = l + 0(z~ l ), z^oo, 

4. T(z) = 0(l), z^O. 

where Jt{ z ) is the following matrix 

/ e M«+W-^_W) e M(e£+(*)-0£_(*)) e M(flX + W-^_W)\ 

J r (z)= e ^(^+W-^-W) 

v e M (^+W-^W)y 

~9f{z)=6f{z)+0UK)- 



(4.11) 



(4.12) 



(4.13) 



By applying Lemma E] to the 8^(z), we can simplify the jump matrix Jt{z). In particular, 
on [A^.AfJ, we have 



Jr(z) 



V 



o 
o 



1 e ^(0£+ (*))\ 

,M(flj v (*)-^(*))_ 



while on [A£_i, A?], we have 



Jt(z) 



\ 










e M(0£ + (*)-«£_M) 
1 






1 



1 





(4.14) 



/ 



\ 



M(eN(z)-0»(z))_ 



(4.15) 



On the rest of the positive real axis, the jump matrix becomes 



Jr(z) 



1 





1 



(4.16) 



This is because McnPn = N% and Mctv(1 — Pn) = N are both integers. And on the 
negative real axis, the matrix T(z) has no jump for the same reason. Note that the jump 
matrix Jt(z) is continuous at z = as the off-diagonal entries of ( 14.16P contain the factor 
e ~M8 1 _(z) ^jgjj vanishes at the origin. 

The Riemann-Hilbert problem for T(z) now takes same form as the one in [TT] (See 
also |27J) and the techniques developed there can now be applied to our problem. 
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Figure 2: The lens contours around a branch cut. 



4.2 Lens opening and approximation of the Riemann-Hilbert 
problem 

We will now apply the lens opening technique that is standard in the Riemann-Hilbert 
analysis. (See e.g. [9], [IE], [17], [18], [TT], j27J) Let us define the lens contours E± around 
a branch cut [A^_ x , A^] as in Figure El We will chose the lens contours such that they are 
contained inside the neighborhood Uj stated in Lemma El Now note that the jump matrix 
Jt{z) in (I4.14p and (14.151) has the following factorizations. 

^(z)J s (z)K J>+ (z), z e [A^_ 1; A^], 




z G [A fe _ 1 , A^J, 



(4.17) 



where Lj(z), Kj(z) are the following matrices 

1 

L 2 (z) = [ eM^W-^M) 1 - e M (^W-^(^)) 

1 

1 

1 

1 ON 
L 3 (z) = | 10 

ve M(0j»(*)-^(*)) _ e Af(^(*)-^W) X 

1 N 

/v. I 1 

,M(0?(z)-9?(z)) e M(§?(z)-9?(z)) x 

and the ± indices in (14. 17ft are the boundary values of the matrices Lj(z) and Kj(z) on 
the branch cuts. 



(4.18) 



28 



If we define the matrix S(z) to be 



T(z), z outside of the lens regions; 

S(z) = { T{z)Kj l {z), z in the upper lens region of [A£J_ l5 A£]; (4.19) 
T(z)Lj(z), z in the lower lens region of [A^_ 1; A^]. 

Then by the factorization (I4.17P and the jump of T(z) (14.161) . we see that the matrix S(z) 
is a solution to the following Riemann-Hilbert problem. 

1. S(z) is analytic in C \ (R + U H J ± ), 

2. S + (z) = S4z)J s (z), ze(M + UHi), 

3. S(z) = I + 0(z- 1 ), z^oo, 

4. S(z) = 0(l), z^O. 

where the matrix Js(z) is given by (14.171) on [X^-i, A^.] and on E±, it is given by 

J s (z)=K 3 (z), zeZ{, J s (z) = L 3 (z), ze-l. (4.21) 

On R + \ (uUK-i, Ag]), we have J s (z) = J T (z). 

Then by Lemma [7] and Lemma [HI we see that, away from [A^ r _ 1 ,A^] and from some 
small neighborhoods Dj of A^, the off-diagonal entries of Js{ z ) are exponentially small 
as M — > oo. This suggests the following approximation to the Riemann-Hilbert problem 

1. S°°(z) is analytic in C \ (u^fA^, Ag]), 

2. = 5?(s)Js(*), ^GU^IAg^.Ag], ( 4 -22) 

3. S°°(2) = J + O^ 1 ), 2^00. 

In the next section we will give an explicit solution to this Riemann-Hilbert problem and 
we will eventually show that S°°(z) is a good approximation of S(z) when z is outside of 
the small neighborhoods Dj of the branch points A^. 



4.3 Outer parametrix 

The construction of the outer parametrix S°°(z) is similar to the one in [TT] (See also [27]). 

Let C N be the Riemann surface defined by (14. 8 h and let Tj be the images of [A^_ 1; A^] 
on £. N under the map £x + (z). That is 

Ii = {(*,£) e£"| £ = *G[A2-i,Ag].}, j = 2,3. (4.23) 
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Let us now define the functions S^(^), k — 1, 2, 3 to be the following functions on £ 



A' 



ST(0 = a A 



a — 1 



1 — a 



nL(£-7j 



n=i(^-7j 



(4.24) 



where 7- are the roots of polynomial 



a 2 (l - c N )^ + 2(a 2 (l - c N (3 N ) + a(l - Cjv (l - #v))£ 3 

+ (1 - Cjv (l - /?at) + a 2 (l - C7V /3tv) + 4a)£ 2 + 2(1 + a)£ + 1 

4 (4.25) 

= a 2 (l-c w )n(e 

3=1 



7f). 



The branch cuts of the square root in ( 14. 24ft are chosen to be the contours Tj (14.231) that 
joins 7^ to 7^ and 7^ to 7^. 

By using the asymptotic behavior of the functions Cf( z ) (13. 13[) and (13.141) . with c and 
P replaced by cjv and (3^, we see that the all the functions S%°(£) are holomorphic near 
£^(0) for j — 1,2 and 3. Moreover, at the points £^(00), these functions satisfy 



^°°(e v (oo)) = ^, kj = 1,2,3. 



(4.26) 



Let us define S°°(z) to be the following matrix on z G C. 



(4.27) 



Then, since the branch cut of the square root in (14.241) are chosen to be Tj, we have, from 
the jump discontinuities of the £f{z) (13.321) . the following 



S?(t?+( z )) = SritZ-i*))' z e [X kj ^, X kj ], 1 = 1,2,3, j = 2,3. 



(4.28) 



From this and the asymptotic behavior ( 14.261) of the 5 , °°(£), we see that the matrix S°°(z) 
satisfies the Riemann-Hilbert problem ( 14.201) . 
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Proposition 4. The matrix S°°(z) defined by \^.21 ) satisfies the following Riemann- 
Hilbert problem. 



1. S°°(z) is analytic inC \ fu| =1 [Ag_ 1; Ag] 

2. S»(«)=5! (z)J5W, ^u^lA^.Ag, 

3. =/ + 0(^- 1 ), 2^00, 

4. O ((*- Af)"*) , J = l,...,4. 



(4.29) 



where Js(z) is defined as in M-llfy . 

Proof. As we have already verified the jumps and the asymptotic behavior at z = oo, we 
only need to verify that properties 1. and 4. in (14.291) are true. From (14.241) . we see 
that as a function in z, the matrix (14.271) has only jump discontinuities along the branch 
cuts [A^_ 1; A^] of ^j v . Also, it can only have singularities at the points 7^ in which the 
denominator vanishes, or at z = in which £i(z) has a pole. As we have already pointed 
out that all the functions Sj°(£) are holomorphic near the point z — 0, the only possible 
singularities are the points z = X^. Near these points, 2 of the functions £f(z) will behave 
as 

= Ik + C s {z - Af)i +0((z- Af )) , * - A?, 

for some constants Cj. Hence the matrix (14.271) will have a fourth root singularity at these 
points. This completes the proof of the proposition. □ 



4.4 Local parametrices near the edge points Aj; 

Near the edge points A^, the approximation of S(z) by S°°(z) failed and we must solve the 
Riemann-Hilbert problem exactly near these points and match the solutions to the outer 
parametrix (14.271) up to an error term of order 0(M _1 ). To be precise, let 5 > and let 
Dk be a disc of radius 5 centered at the point A^, k = 1, . . . , 4. We would like to construct 
local parametrices S k (z) in Dk such that 

1. S k (z) is analytic in D k \ (R U E k _ U S*) , 

2. Si (z) = S k _{z)J s {z), z G D k n (R U S fc U S*) , (4.30) 

3. S k (z) = (/ + 0(M- 1 )) S°°{z), zedD k . 

The local parametrices S k (z) can be constructed by using the Airy function as in pT] (See 
also [27]). Since the construction is identical to that in [11] and [27J, we shall not go into 
the details but merely set up the notations and state the results in [TT] and [2"T] . 

First recall that the Airy function Ai(z) is the unique solution to the equation v" = zv 
with the asymptotic behavior given by (11.191) . 
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Let uj = e z\ then functions Ai(z), Ai(uz) and Ai(u 2 z) satisfy the following linear 
relation, 

Ai(z) + loM(ujz) + Lu 2 Ai(u 2 z) = 0. (4.31) 
Now note that, in the neighborhoods D k and D k -1, j = 2, 3, the functions 



f*(z) = e?(z) - vCAg) - ef{z) + %,_(Ag), 



(4.32) 



3 , .3 



vanish like — Aj^j 2 and ^2 — ^kj-ij 2 respectively, as z approaches or A^_ x from 

the lower half plane. Let us take the restriction of ff'iz) and ff{z) to the lower half plane, 
and analytically continue them to D k . \ [X%. , A^ f _ 1 ] or to D^—i \ [X k . , A^_ x ] . We will denote 
these analytic continuation by f^{z) and fh(z). Then the functions fjj(z) and fh{z) will 

behave like — A^j 2 and (^z — \ k ,_^ 2 in D k . and Dk 3 -i respectively. Therefore if we 



define the local function ((z) by 



2 

4 \ s 



3 ' (4-33) 

then £(z) will be holomorphic inside the neighborhoods D k . Moreover, for small enough 
5 > 0, ((z) will be conformal inside D k . (The | factor is introduced to simplify calculation) 
We should make use of the freedom in the definition of the lens contours to deform them 
such that inside D2 and D4, the upper and lower lens contours coincide with the contours 
axg(((z)) = ^ and arg(£(z)) = — respectively, while inside D\ and D 3 , we will deform 
the lens contours such that the upper and lower lens contours coincide with arg(£(z)) = | 
and &rg(C(z)) = — f respectively. 

Let us now consider the matrix S (z) that is related to S k (z) by 

z inside the lens of Dk 2 -u 
S k2l (z)=^ S k 2 -i (z) ( 1 _ e M(e»(z)-S»(z))\ z outside the lens of L^. 



S^-Hz) 




z inside the lens of D^-u 



(4.34) 




z outside the lens of D 



ki - 



where i = 0, 1 in the above. In [TT] (See also [27J), matrices S k (z) was constructed such 
that when S k (z) is related to S k (z) by (Q3|) . S k (z) will solve the problem (IPO]) . These 
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matrices have the following form (See also [9], [IE], jTTj , [T8]). 

S k (z) = E k (z)W k $ k (m§C(*)) ^ G 
where Z k is the diagonal matrix 

Z k = diag (eiWW), c -i"(*f M-arM), i) , = fc 2 , A; 2 - 1, 



(4.35) 



Z fc = diag (e^^W-^W), 1, g-i^W-^W)) , k = k 3 , k 3 



1. 



(4.36) 



The matrix E k (z) is a holomorphic matrix inside that is bounded as M, N and N\ — ► oo. 
The matrix is a constant invertible diagonal matrix of order Me. They are here to 
fix the boundary condition in (14.301) . The Q k (z) is a matrix whose entries consist of the 
different branches of the Airy function. 

For example, when k = 4 and a < 1, the matrices E k (z), W k and & k (z) are given by 



%(z) -v 2 (z) N 
v' (z) -v' 2 (z) 
1, 

-«i(z) -v 2 (z) N 


1, 

-v 2 (z) ^(z) N 
-t£(z) v[(z) 
1, 
%(z) v x {z) N 

1 



for < arg(z) < 



2tt. 

3 ' 



for =^ < arg(^) < 7r; 



for —it < arg(z) < 



for 



— < arg(z) < 0. 



EAz) 



Wa 




diag(M^,M~^,l) 



where the functions Vj(z) are given by the Airy functions Vq(z) = Ai(z), V\[z) = ujA\{ujz) 
and v%(z) = u 2 Ai(uj 2 z). Note that despite the apparent singularity in the expression of 
E±(z), it turns out that E 4 (z) is holomorphic inside D± because the factor S°°(z) contains 
terms that would cancel out this singularity. 



4.5 Last transformation of the Riemann-Hilbert problem 

Let us now show that the parametrices we constructed in Section 14.31 and Section 14.41 are 
indeed good approximation to the solution S(z) of the Riemann-Hilbert problem (14.20]) . 



33 



Figure 3: The contour Tr. 



Let us define R(z) to be the following matrix 

r,f z) = \ S{z){S k {z)Y l , z inside D k , k = 1, . . .,4; , , 

y> \ S(z)(S O0 (z))- 1 , 2 outside of D k , fc = l,...,4. l ' J 

Then the function -R(^) has jump discontinuities on the contour Tr shown in Figure El In 
particular, R(z) satisfies the Riemann-Hilbert problem 

1. R{z) is analytic in C \ Tr 

2. R + (z) = R-(z)J R (z) 

3. R(z) = I + 0(z- 1 ), z^oo, 1 ' ' 

4. R(z) = 0(1), z -> 0. 

From the definition of (14.371) . it is easy to see that the jumps Jr(z) has the following 
order of magnitude. 

j (z) = f 1 + OiM- 1 ), zedD k ,k = i,...A; f4 39) 

R ^ ' \ I + O (e~ Mr >) , for some fixed 77 > on the rest of T R . k ■ J 



R{z) =1 + 01 —rr-^, ) , (4.40) 



Then by the standard theory, [IB] . [T7] . [18], we have 

% M(jz] + 1 
uniformly in C. 

In particular, the solution S(z) of the Riemann-Hilbert problem f)4.2Up can be approx- 
imated by S°°(z) and S h (z) as 

q( 7 \ = \ (I + O (M^ 1 )) S k (z), z G D k , fc = 1, . . . , 4; 
l J \ (J + 0(M- 1 ))5 0O (2), 2 outside of D k , k = 1, . . . , 4. 1 } 

5 Universality of the correlation kernel 

We shall now use the asymptotics of the multiple Laguerre polynomials obtained in the 
last section to prove the universality results in Theorem [3j Since the proofs are the same 
as the ones in Section 9 of [TT], (See also [9], [16], [TF], [H]) we shall only carry out the the 
proof for (11. IT)) and leave the readers to verify (jl.isp . 
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First let us express the correlation kernel Km,n{x,u) in (11.151) in terms of the matrix 
S(x) in 105)1 . By g2J), ffTToD . fTCTTTl and (ICTjl . we have 

K MA^v)= ^_ 2 y) (-e Me ?. + M e M§ ^ 0)S-\y)S + (x) L-^Wj , (5.1) 

for x,y (z [Ag_!,Ag] and 

_ L j^iv / e -MO» + (x)\ 

K M A^y)= [ X ' ( y) \ U M ^ e A ^)) S;%)S + (x) (5.2) 

for x,y (z [A^_ 1; A^]. Now note that, since Ci( z ) an d ^fi 2 ) are complex conjugates on 
[Ag_i,Ag],byi3J,wehave 



< + W-< + (A^) = ^ + W, (5.3) 



where A denotes the complex conjugation. 

If kj = 4, then 0f + (A£) = 0i,-(A£) = and if fc, = 2, then 



^±(AD = / ZU*)** 



A? 
Af 



Since = £%_(z) in [Af ,Af] and = £%_(z) G R in [Af , Af ], we see that, in 

either case of fcj = 2 or fcj = 4, we have 0f + (Aj^) = 0f_(Af ). Hence we have, by (15.31) and 

USD, 

= 0£ + {z), ^[A^AQ. (5.4) 

By the same argument, same relations between Oj(z) and 9j(z) can be obtained. 
By substituting (15 .4p and (I1.15P into (15.11) . we obtain 



2iri(x — y) 

/ e -MIm(e» + {x))\ 

V / 

on [Aj^_ x , A^], where hjsr{x) is given by 



ii MAx , y) = e-""^>-"" W > e -«-« + M) o) «:'(») 



x S+(x) 



(5.5) 



M*0 = Re « + (x)) - ~ (0 li+ (x) + 9 h+ {x) + 0i,_(A fc .)) 
= Re(^ + (x)-0 1 , + (x)). 

35 



where we have used ( 15.41) in the second equality. We can obtain similar expression for 
x,ye [Aj_ ls Ag] by substituting (E3D into (E2D. 

We can now apply (14.411) to (15.51) and its counterpart on [A^_i, X^] to prove (11. 17ft in 
Theorem [3J 

Let us consider the points x and y G (A^-i, A&.) such that 

u v 
Mp N (x ) Mp N (x ) 

where xo G (A^._i, A^.). For large enough N, we can assume that xo, x and y are in 
(A^_i, Aj^) and outside of Dk 3 -i and D^.. The function Pn(z) in (]5.6p is given by 

Pn{z) = -Im , (5.7) 

Let us replace the matrix S+ (y)S + (x) in (15.51) by 

S-\y)S + (x) = {S^iy))- 1 R- 1 (y)R(x)S^(x). 

Since -R(z) is analytic at xq, by (I4.40H and by considering the power series expansion of 
R(z) at xo, we see that i? _1 (y)-R(x) is of order I + (^jf)- On the other hand, since xq 
is outside of the neighborhoods of the branch points Aj^, the matrix- valued function 
ST(z) has an analytic continuation in a neighborhood of Xq that includes both x and y. 
Hence by considering the power series expansion of S™(z) at xq, we have 



S- l (y)S + (x) = (sriy))' 1 R-\y)R(x)S^(x) 



— I + 0(x — y) 

By substituting (15. 8p back into (15.51) and its counter part on [A^_ 1 , A^], we obtain 

K MjN (x, y) = e -M(h N ( x )-h N{y) ) V 1 ' +l \ h+yy))J + 0(1) , (5.9) 

\ n(x-y) J 

for 1,1/6 (A/^_i, A^). Since 9^ + (z) is given by (14.91) and that xo is away from the branch 
points A^, the functions 6± + {z) and £i + (z) can be continued analytically in a neighborhood 
of xq. If we denote these analytical continuation by the same symbols, we will have 

= CfM + Ci,+{ z ~ x °) + (i z ~ x of) , z^z . (5.10) 
Now by substituting (15.61) and fl5.1Up into (15.91) , we obtain 



sin i Im I £iidil_!l2 | — v)\ 
K M , N (x,y) = e -"<M»)-MW) | MpN (x ) V V P » M '1 + 0(1) ) . (5.11) 
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From the definition of Pn(z) 115. 7ft . we see that 

Im — — — - | 7T 

\ Pn{x ) 

Now recall that the density p(z) on the support of F_ is given by (See Theorem [U]) 

p(z) = ilm (&,+(*))■ 

7T 

As the solutions of ( 13. 12ft are continuous in the parameters a, (3 and c, we have 

lim ^(zo) = £i,+ (x ), Jim piv(x ) = p(x ). (5.12) 

Let us now show that 

lim e -mN(x)-h N (y)) = L / 5>13) 

M,N^oo 

By dSSD and (15~TUD . we see that 

lim M (ftjv(ar) - h N (y)) = lim Re (^+(x ) - Zi,+ (x )) — 
Then by CTl . we obtain (15131) . 

1 1 h Vttr 4~Vin toqI itq ri q l~iloe 11 — i I ( 



Now if we replace u and v in (I5.1ip by the real variables u — > " p y^° , u — > vp ^ x °^ and 



take the limit, we obtain 



sin 7T [u — V) 



lim - — = p(x ) z — L . (5.14) 

This proves (11.171) . 

The proof for ( 11.181) is the same as the one in Section 9 of [TT]. Let x and y be 

x = Xk+ r, l/ = Af + r . (5.15) 

(p»M)3 {p N k MY 

Then following through the details in Section 9 of [11] . we would arrive at 



-M(h N (x)-h N (y)) / ^ 

Km,n \ H —^2i A k 



(Mp»y V (Mpf)3 {MflYJ ( 5.i 6) 

Ai(u)Ai» - Ai'(u)Ai(i;) 



u — V 



0[M~* 



We will leave it to the readers to verify the details. To obtain ( II. 18ft from ( 15. 16ft . we will 
need to use the following. 
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Lemma 9. Under the limit A& satisfies the following 

\ k = \% + (M" 1 ) 

Proof. Let a e be a solution of the following quartic equation 

W (a e ) + eW 1 (a e ) = 0. 

where Wj(x) are polynomials of degree n — j whose coefficients are bounded in e. Let «o 
be the solution when e = 0. Suppose ot e = a + 5. Then we have 

(Wo(a ) + eWi(a ) + 0{5 2 ) + 0(Se) + 0(e 2 ) = 0. 



Therefore if | = o(l), we will have 



Wo(a ) = 0. 



This implies that ao is a double root of Wq(x). We can apply this to the polynomial (11.8)1 
and (T05D . By (Oil , we have c - c N = O (M -1 ) and /3 - /3jy = O (M -1 ). Therefore the 
differences between ( II. 8p and (I4.25P is a cubic polynomial whose coefficients are of order 
O (M^ 1 ). 

Since the polynomial ( 11.81) has no double root, the difference between A& and X k must 
be of order O (M _1 ). □ 

Let us now show that ( 15. 13ft is true in this case. Near the edge point A&, the function 
h N (x) - h N {y) is given by 

3 

2 I U2 — V2 \ 2 I I , at , u 



h N (x) - h N (y) = - ^ — — — j - -p k ( ( A A - A/, -! 
2 



3 VV (p£M)*JJ v 7 

Since A^ — A^ = O (M _1 ), the above becomes 

M.)-Mv) = f(i-£) (%^) +0 (m-* 

This implies that (15. 13[) is true also at the edge. 
Now let us replace u and v in (15.161) by 

u-(^)^-(^M)MAf-A & ), 



^ «-(rfM)»(Af-A fc ). 

Since Aj^ — X k = O (M _1 ), the new u and v are still real and of the same order in M. By 
making these replacement, and then take the limit of (I5.16p . we arrive at (11. 18j) . 
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